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ABSTRACT 

Weak lensing by large-scale mass inhomogeneities in the Universe induces correlations 
in the observed ellipticities of distant sources. We first review the harmonic analysis 
and statistics required of these correlations and discuss calculations for the predicted 
signal. We consider the ellipticity correlation function, the mean-square ellipticity, 
the ellipticity power spectrum, and a global maximum-likelihood analysis to isolate 
a weak-lensing signal from the data. Estimates for the sensitivity of a survey of a 
given area, surface density, and mean intrinsic source ellipticity are presented. We 
then apply our results to the FIRST radio-source survey. We predict an rms ellipticity 
of roughly 0.011 in 1° x 1° pixels and 0.018 in 20' x 20' pixels if the power spectrum 
is normalized to as^'^'^^ = 0.6, as indicated by the cluster abundance. The signal is 
significantly larger in some models if the power spectrum is normalized instead to the 
COBE anisotropy. The uncertainty in the predictions from imprecise knowledge of the 
FIRST redshift distribution is about 25% in the rms ellipticity. We show that FIRST 
should be able to make a statistically significant detection of a weak-lensing signal for 
cluster-abundance-normalized power spectra. 
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1 INTRODUCTION 

It has been proposed that the effects of weak lensing on dis- 
tant sources could shed l ight on the l arge-scale distribution 



Blandford et al. 199 1[ Kaiser 1992|; 


3artelmann & Schneidei 


1992; Kaiser 1996; Villumsen 1996 


; Stebbins 1997; Jain & 



Seljak 1997). Mass inhomogeneities along the line of sight 
to distiiiiL Huurtiya will Indutiu dlyLonioiiy In Lliy IniagoH of 
these souicys, thus, LuiiylaLiuuh ul Llie ellip Liu Lies uf disLauL 
sources provides a probe of the correlation of mass along 
the Une of sight. In this way, the power spectrum for the 
mass (rather than light) distribution in the Universe can be 
probed. 

This technique requires (i) a sample of sources which are 
distant, so there is a large line of sight over which the lensing 
signal can accrue; (ii) good angular resolution, so that the 
ellipticities of sources can be determined; and (iii) a large 
enough sample of sources so that the noise provided by the 
intrinsic ellipticities of sources can be overcome. For exam- 
ple, the VLA FIRST radio survey m eets all of these criteria 
( Becker et al. 1995 ; White et al. 1997 ). Upon completion, the 



survey will cover 10,000 square degrees of the North Galactic 
cap. There are ~ 40 sources per square degree with resolved 
structure on scales from 2 — 30" at the survey resolution 
of 5", and the mean redshift of these sources is of order 
unity. Several systematic effects can produce spurious ellip- 
ticity correlations in FIRST, and therefore be mistaken for 
a weak-lensing signal. However, the most seriou s systematic 
effects are understood and can be corrected for (Refregier & 



Brown 1997 



Refregier et al. 1997) 



When looking for the effects of weak lensing on galaxies 
behind a cluster, one requires tens of thousands of galaxies 
per square arcminute to overcome the Poisson noise from in- 
trinsic galactic ellipticities and thereby map the shear field. 
With this in mind, it may at first seem hopeless to detect the 
effects of weak lensing in a sparsely-sampled survey such as 
FIRST with < 100 resolved sources per square degree. How- 
ever, for weak lensing from large-scale structure, we are in- 
terested in the correlation of ellipticities of pairs of sources 
with some fixed angular separation; we do not necessarily 
need to map the shear field. For this, the relevant quan- 
tity is not the density of sources, but the total number of 
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pairs of sources in the survey with some fixed angular sep- 
aration. In other words, the sensitivity to the mean-square 
ellipticity averaged over regions of some fixed size on the sky 
is improved with a large area [cf., the discussion of sparse 
sampling in Kaiser (1996)]. 

In this paper, we review the theory of ellipticity correla- 
tions from weak gravitational lensing. We discuss statistical 
t echniq ues which can be used to isolate a signal in the data. 



We also estimate the amplitude which may be detectable 



with a survey which covers a given fraction of the sky with 
a given number of resolved sources per square degree and 
a given mean intrinsic source ellipticity. We then calculate 
the predicted weak-lensing correlations in the FIRST radio 
survey for the canonical cold-dark-matter model as well as 
for several viable variants. We discuss the uncertainties in 
the predictions which arise from imprecise knowledge of the 
redshift distribution of FIRST sources, and consider the de- 
tectability of the signal. Refregier & Brown (1997) show how 
spatial noise correlation, one of the most serious systematic 
effects for FIRST, affects weak-lensing measurements and 
can be corrected for. A subsequent paper will pr esent the 
results of the search for weak lensing in FIRST (Refregier 
et al. I!p97|). 



2 ELLIPTICITY CORRELATIONS FROM 
WEAK LENSING 

2.1 Description of the Shear Field 

Weak lensing will induce a stretching of images on the sky 
at position d = (dx, 9y) described by the shear field, a sym- 
metric, trace-free 2x2 tensor field: 



(1) 



Here, e+ is the stretching in the 9x — Qy directions, and ex is 
the stretching along axes rotated by 45°. Alternatively, the 
shear field can be written as a shear "vector," 



7q = (e+,ex) = 7(cos2a,sin2Q), 



(2) 



although this ordered pair does not transform as a vector. 
The deformation is also sometimes written as a complex 
ellipticity p; the two components of the polarization "vector" 
are the real and imaginary parts of the complex ellipticity. 



P = 



\p\e 



(a2 + 62) 

/ 2 I 2 \l/2 2ia , • 



(3) 



where a and b are the principle moments. We see that 7 = 
IpI = (el+el)^^^ = {a^-b^)/ia'^ + b'^) and tan 2a = e^/e+. 

As pointed out by Stebbins (1997), the shear tensor field 
laffiO), can be decomposed into a "gradient" or curl- free 
component (what Stebbins calls the scalar part) and a "curl" 
(or pseudoscalar) component, just as a two-dimensional vec- 
tor field can be decomposed into curl and curl-free parts. In 
other words, the shear tensor can be written in terms of two 
scalar functions 79(6') and 7c (6*), 

(4) 



where ea/} is the antisymmetric tensor. Specification of 7g (6*) 
and 7c (^) is equivalent to specification of £+(6*) and 



The quantities 7g and 7c are the gradient and curl compo- 
nents of the ellipticity field. 

Density perturbations (mass inhomogeneities) produce 
only scalar perturbations to the spacetime metric, so they 
have no handedness and can therefore produce no curl. 
Gravitational waves have a handedness and may induce a 
nonzero 7c. However, the weak-lensing signal from gravita- 
tional waves is expected to be extremely small ( ^tebbim 
1997). Vector modes could also produce a curl, but, like 



tensor modes, they are negligible. Therefore, we expect that 
7c = 0, and only 7g should be nonzero. This allows one to 
isolate the effect of weak lensing and to check for non-lensing 
ellipticity correlations. 

Throughout, we approximate the region of sky surveyed 
as a flat surface. At first this might seem inappropriate since 
the FIRST survey covers a good fraction of the sky. However, 
weak-lensing ellipticity correlations should be most signifi- 
cant only at smaller angular separations, so the inaccuracies 
which arise from approximating the sky as a flat surface 
should be small. Furthermore, we are primarily interested 
here in establishing the existence of an ellipticity correla- 
tion. With future more sensitive data, i t will be neces sary 
to account for the curvature of the sky (Stebbins 1997). 

Given the shear tensor 7^/3 (S), the functions 79 and 7c 
can be constructed only with a Fourier transform. Writing 

j2,- 



(2^ 



7«/3(k)i 



(5) 



where the tilde denotes the Fourier transform (and similarly 
for the other quantities), we get 



7s (^) 



7c (k) 



{ni — kI)£+{k) + 2KxKyex (k) 



(6) 



(7) 



The functions 7g(6) and 7c (6) can then be recovered through 
the inverse Fourier transformation. Since 7^^ is a real tensor, 
7q/3(^) = 7(1/3 (—k), and similarly for 79 and >. 



2.2 Power Spectra 

Statistical homogeneity and isotropy guarantee that the two 
sets of Fourier coefficients, 79 and 7c, have expectation val- 
ues. 



(79(^)79('^')) 

(7c(^)7c(k')) 
(79(^)7c(k')) 



(8) 



The power spectrum Pgg{K) is precisely Kaiser's (1992) el- 
lipticity power spectrum Pe{K,). The second power spectrum 
Pcc(k) will be effectively zero because weak lensing from 
gravitational waves is extremely small. The third, Pgc{Mi), 
must be identically zero since it is parity violating — that is, 
this power spectrum changes sign under the change of co- 
ordinates X —> —X. Since these latter two power spectra are 
zero, they can be used to look for non-lensing artifacts in 
the data. 
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The mean-square gradient component of the eUipticity 



(27r)2 



(9) 



Since the mean-square curl component of the eUipticity is 
zero, (79) is also the mean-square total eUipticity. 

Realistically, the mean-square eUipticity cannot be mea- 
sured. The actual measured quantity is the mean-square el- 
lipticity smoothed with some window function. Suppose we 
estimate the shear field at position 9 by averaging over all 
ellipticities, e.g., in a square Op x Op pixel centered at 0. In 
that case, we are probing a smoothed shear field. 



4W 



(10) 



where W{a) is the window function (e.g., for square pix- 
els, constant inside the pixel, zero outside, and normalized 
to unity) , and similarly for e ^ • prior results for the un- 
smoothed field are generalized to the smoothed field as long 
as we replace £+(«:) 2+('^) = ^+(^)W(Jl), where Win) is 
the Fourier transform of W{9). The mean-square smoothed 
eUipticity is then 



(11) 



In the following, we also use the shorthand Pgg(hi) = 
Pgg{K)\W{K)\^ for the smoothed power spectrum. Although 
the mean-square eUipticity gives a simple indication of the 
magnitude of the weak-lensing signal, one can obtain a much 
more sensitive probe of a signal by taking advantage of the 
information provided by the complete power spectra (or 
equivalently, correlation functions), as discussed further be- 
low. 



2.3 Correlation Functions 

There are three independent two-point eUipticity correla- 
tion functions that can be constructed in configuration space 
from the three power spectra. Since the components (e+ and 
ex) of the shear tensor are defined with respect to some set 
of axes on the sky and transform under rotations, correlation 
functions of these quantities will depend on the relative ori- 
entation of the two points being correlated as well as the sep- 
aration. However, correlation functions which are indepen- 
dent of the coordinate system can be constructed ( Stebbins 



1997) in ana logy with those needed for CMB polarization 
correlations ( Kamionkowski et al. 1997). To do so, we de- 
fine correlation functions of ellipticities e!^ and measured 
with respect to axes which are parallel and perpendicular 
to the line connecting the two points being correlated. To 
be explicit, suppose the first point is Oi — {0ix,Oiy) and the 
second is 02 = {02x,O2y). Then we must rotate the axes by 
an angle = arctan[(6'2y — 0iy)/{02x — Oix)] to align the ro- 
tated X axis with the line connecting the two points. Under 
this rotation, we get 



e_i. = cos 20 + Ex sin 20, 
e'x = — e+ sin 20 -f ex cos 20, 



(12) 
(13) 



at both points. The 20 enters since the eUipticity is un- 
changed under a rotation by 90°. We can then construct 



three correlation functions, {e!j.e!^), (e'x ex), {elj.ex), from the 
rotated components. Although e+ is invariant under reflec- 
tion along the line connecting the two points being corre- 
lated, ex changes sign. Therefore, parity invariance demands 
that (e+ex) ~ 0. Statistically significant deviations from 
zero can be due only to systematic errors in the data. 

By setting = in Kaiser's Eq. (2.3.1), we iden- 
tify (e;(eo)e!;(^o + e"')) = Ci(e) and (e^AOoyxido + 0)) = 
C2{0). We also verify that (e!;(eo)e^ {0o + 0)) = 0. In anal- 
ogy with correlation functions of Stokes parameters of th e 
cosmic microwave background (Kamionkowski et al. 1997), 
we can write the correlation functions (for any 0) in terms 
of the power spectra as 



Ci{0) + C2{0) = 
a{0)-C2{0) = 



^[P,,(A.) + P..(^)]Jo(«e) (14) 



Pcc{>^)]Ji{K0). (15) 



In terms of correlation functions of rotated and unrotated 
ellipticities, and in terms of the complex eUipticity p, 

C,(0) + C2{0) = {e\{0o)e\{0o + 0)) + {eAeiMx0o + 0)) 
= {e+{0o)e+{0o + 0)) + (ex {0o)ex {0o + 0)) 
= Re[p*{0o)p{0o + ff)], (16) 

Ci [0) - C2{0) = (e; {0o)e\ (00 + 0)) - (e; {0o)el {0o + 0)) 



cos 40 



{e+i0o)e+{0o + 0)) 



(ex(eo)ex(^o + 0)> 



4- sin 40 



(ex(^o)e+(^o + e)) 



+ (e+{0o)e^{0o + 0)) 



Re[p(6'o)p(6'o + e)]cos40 
+lm[p{0o)p{0o + 0)]smAq 



(17) 



where — cos 0, 6* sin 0}. There is also the third linearly 
independent correlation function, 

CcrosB [0) = ( e'+ (^0 ) e; (^0 +0)) + { 4 (^"o ) (^"o + 0) ) 

= - sin 40 [{e+{0o)e+{0o + 0)) - (ex {0o)e^ {0o + t 
+ cos 40 [(ex (^o)e+ (^o + ^)) + (e+ (^o)ex (^o + ( 

= -Re[p{0o)pi0o + 0)] sin 40 
+Im[p(^o)p(9o + e)]cos40, 

and parity conservation demands Ccross(^) = 0. Note that 
when written in terms of the wnrotated ellipticities or the 
complex eUipticity, the sum Ci{0) +C2{0) is independent of 
0. However, when written in terms of the unrotated elliptic- 
ities or complex eUipticity, the difference Ci{0) — C2{0) and 
Ccross(6') does depend explicitly on 0, the relative orientation 
of the two points being correlated. 

The power spectra can be written in terms of the cor- 
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relation functions as 



+[Ci{e) - C2{e)]MKe) 



-[Ci{e) - C2(e)]MKe) 



(19) 



Again, if the second of these is nonzero, it can only be due to 
non-lensing effects, so construction of this correlation func- 
tion provides a powerful probe for the presence of nonlensing 
artifacts in the data. 

Of course, correlation functions C (d) for the smoothed 
ellipticities can be obtained by replacing ellipticities and 
power spectra by the smoothed ellipticities and power spec- 
tra in all the equations above. 



3 STATISTICAL ESTIMATORS 

The effects of weak lensing can be uncovered through the 
measured correlation functions, power spectra, mean-square 
ellipticities averaged over some given pixel size, or a full 
maximum-likelihood fit to the data. 

The Fourier modes of the shear field due to weak lens- 
ing are statistically independent. Furthermore, if the noise 
map is orientation independent, then the Fourier modes of 
the noise will also be statistically independent. Even if we 
use only the simplest (although not necessarily optimal) es- 
timator for the power spectrum, the mean-square ellipticity, 
it is better to work in Fourier space. The predicted weak- 
lensing mean-square ellipticity is due entirely to the gradient 
component, but randomly oriented intrinsic source elliptici- 
ties should contribute to the mean-square ellipticity equally 
through the gradient and curl component. The signal-to- 
noise ratio will therefore be improved with a Fourier trans- 
form which allows us to isolate the gradient and curl com- 
ponents. 



3.1 Discrete Fourier Transforms and Statistical 
Noise 

We restrict ourselves to a survey which covers a rectan- 
gular region of the sky. The analysis can be extended to 
irregularly-shaped regions of the sky, but only with signif- 
icant complications. [Simple estimates of the effects of an 
irregularly-shaped survey w hich are used for the power s pec- 
trum of angular clustering ( Baugh fc Efsthathiou 1994 ) are 
not easily generalized to weak-lensing power spectra.] We 
first construct pixels of size Op x 6p on the sky. This leaves 



us with Np 



Nx X Ny pixels, where — 6., 



■JOp 



and 



0y/9p and O^. and 9y are the dimensions of the map. 
The pixels labeled by (i, j) are centered at 9ij — {i,j)8p, and 



0,l,...,iV,- land; =0,1, 
The ellipticities, e 



,Ny 



1. 

measured in pixel 



are the sum of a weak-lensing signal and noise which arises 
from intrinsic source ellipticities and measurement error. 



+ e 



(20) 



We can then use an FFT to determine the A'^pix Fourier 
coefficients of the survey. 



~obs / -*\ X ^ obs ii 

ij 



-obs / -*\ X ^ obs i 



for fi = {2TT/ep){n/Nx,m/Ny), and n = 0, 1, 



(21) 



1 and 



m = 0, 1, ...,Ny - 1. 



If the noise terms are all statistically independent 
with variances [i.e., they satisfy {^+^ij^^,ki) = (^eSikSji, 

{<^l,ij^Z,ki) = and (e+,yex,fci) = 0], then estima- 

tors for the mean-square ellipticities are given by 



(71)^ 



(7l)*7l 




p;c(«). 



(22) 
(23) 
(24) 



These are estimators for variances of a distribution measured 
with a finite number of pixels. Therefore, there will be some 
cosmic variance as well as some pixel-noise variance with 
which these estimators will recover their expectation values. 
These variances are 



2 



[((7^)^)+a?]^ (25) 



(26) 



iVo 



X [((7c)')+^?] . (27) 

These results may be obtained in analogy with the deriva- 
tion for cosmic and pixel-noise variances for a temperature- 
polarizatiori map of the cosmic microwave background 
( tKnox 199E| ). 

Inserting the null hypothesis of no signal, ((7g)^) = 
0, into Eq. (kaj gives us the statistical limit to the weak- 
lensing amplitude of this quantity to which we are sensitive. 
Explicitly, we can be assured a 3a detection of (7^)^ only if 
it exceeds 3(Jj y'2/iVpi^. 

If the density of resolved sources on the sky is n (in units 
of deg~^) and the mean intrinsic ellipticity of the sources 
is e (the mean intrinsic value of |p|), then = €'^/{n6p). 
Therefore, if the area of the survey is A, then the survey 
will be sensitive (at la) to a mean-square ellipticity. 



-1/2 



'^{(7|)^) = (0.0075)" (^/lO, 000 deg 

x(e/0.4)^(n/40deg"^)"^(6lp/deg 



(28) 



for pixels of area dp. Since the signal is the mean-square el- 
lipticity (rather than the rms ellipticity), an rms ellipticity 
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> 73(0.0075) ~ 0.013 in 1° square pixels should be de- 
tectable at 3(T with the survey parameters assumed here. 
The central values above were chosen to be close to those 
expected for resolved sources in the completed FIRST sur- 
vey. 



3.2 Likelihood Analysis 

Although the mean-square ellipticity per pixel provides a 
simple estimate of the sensitivity of a given survey to a 
signal, it is not the optimal statistic for detecting a weak- 
lensing signal. The sensitivity of a survey to a weak-lensing 
signal can be improved significantly with a maximum- 
likelihood analysis which compares the complete power spec- 
trum (rather than just the mean-square ellipticity) with the 
entire survey. 

Suppose that in a survey with A^pix pixels we 
construct a 2A'pix-dimensional data vector, D"^^ = 

r obs obs obs obs obs obs 1 frnm +>io 9 ^/' • 



measured ellipticities and e+.ij, and each observed 

value is due to signal and noise, D"^^ — + D^. Sup- 
pose further that we are testing a Gaussian theory which 
predicts expectation values (^D^D^^^ = C^p (where the the- 
ory correlation matrix is given by the unrotated correlation 
functions discussed in section 2.3) with a map which has 
Gaussian noise with a correlation matrix {^D^D^'j = C^^. 
The likelihood of this theory given the data is 



rcxexplD^^KC'+C") 



(29) 



For example, if the noise is due entirely to intrinsic source 
ellipticities, then the noise in each pixel is uncorrelated and 
the noise between -|- and x ellipticities is also uncorrelated, 
so the noise correlation matrix becomes diagonal with en- 
tries equal to the variance in each ellipticity, = a^Sap. 
In general, however, the noise correlation matrix will be non- 
diagonal (Refregier & Brown 1997), and the theory matrix 
is also nondiagonal. Therefore, for a 10,000-deg^ survey with 
20' X 20' pixels, the data vector will have 180,000 entries, 
and evaluation of the likelihood would require inversion of a 
180, 000 X 180, 000 matrix! 

Progress in evaluating the likelihood with good accu- 
racy can be made by working in the Fourier domain in- 
stead. In this case, we write the data as a 2A^pix-dimensional 
vector with the A'pix measured Fourier components 7g'"'(«:) 
and 7°'"'(«:) as components. Statistical isotropy and ho- 
mogeneity guarantee that these have expectation values 

(7|(«)7s(^')) = PsA^)Sax'^ (7|(^)7c(^')) = 0, and 
(7°(k)7^(k')) — 0. In other words, the covariance of the 
Fourier components is diagonal, and since all the c compo- 
nents are expected to be zero, we only need to deal with 
a A'pix X Apix (rather than 2Apix x 2Apix) covariance ma- 
trix. Unfortunately, the covariance of the noise in Fourier 
space will not always be diagonal, unless the noise correla- 
tion functions depend only on the separation between two 
points and not their relative orientation. This will not neces- 
sarily be the case. For example, the noise in interferometric 
images is generally spatially correlated and has preferred 
directions (visible as "stripes" in the noise). This produces 
spurious ellipticity correlations which d epend on the orien- 



tatio n of the source-separation vector (Refregier & Brown 
1997). Therefore, the error obtained by assuming the co- 



variance in Fourier space to be diagonal will underestimate 
the true error. By carrying out a full likelihood analysis on 
a small patch of the survey and comparing it with the re- 
sults of the restricted analysis (i.e., that which assumes the 
Fourier components are all statistically independent) on that 
same patch, one can determine the degree to which the re- 
stricted analysis underestimates the error. 

To proceed, we must make the approximation that 
the covariance of the noise is diagonal in the Fourier do- 
main. Then the power spectrum P^g{n) of the noise can be 
obtained from the noise correlation functions, Ci{9) and 
C?(^)i which are obtained by averaging over all orienta- 
tions (j>- With this approximation, both the noise and the- 
ory covariance matrices are diagonal in the Fourier domain, 
and evaluation of the likelihood reduces to evaluation of the 
usual x^- Therefore, if our model consists of a smoothed 
power spectrum APgg{K,) and we are trying to fit for the 
amplitude A, then each of the Apix measured Fourier com- 
ponents K provides an independent estimator. 



Ak 



^(^)|^-P;,(A.) 



for A, with a variance 



2 



, = {{A,-{A,}f) = 



Therefore, an estimator for A for the entire survey is 



A = 



and the variance to this estimator is given by 



1 



(30) 



(31) 



(32) 



(33) 



To determine the smallest value of A which could be 
distinguished from the noise, we insert the null hypothesis 
^ = into Eq. (^). To illustrate, suppose that the noise in 
each component of the ellipticity was Gaussian with variance 
and independent of each other component, as assumed 
above. In this case, Pgg{K) = cr^. Therefore, 



i = 2^Ei^«v^)]^ 



A^. 

' pix 

2af 



(27r)2 



or in other words 



1A. 
A 

where 



1/2 



Ap 



fiA ' 



-1/2 



A' 



[P^,i^) 



(34) 



(35) 



(36) 



(27r)2' 

and we have used — /{ndp) and 6^ = A/Npi^. Eq. ( |35[ ) 
shows how the signal-to-noise scales with the mean intrinsic 
source ellipticity e, usable density of sources n and the survey 
area A. 



4 PREDICTIONS FOR THE FIRST SURVEY 
4.1 Preliminaries 

We will restrict our analysis to a flat Universe (fio + f^A = 1), 
but will allow for a nonzero cosmological constant SIa. The 
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scale factor of the Universe, a{t), satisfies the Friedmann 
equations, 



- = HoE{z) = i/oV^^o(l + z)3+nA, 



(37) 



where Ho ~ 100 /i kmsec"^ Mpc~^ is the Hubble constant, 
Qo is the current nonrelativistic-matter density in units of 
the critical density, Qa is the current contribution of the 
cosmological constant to closure density, and the dot denotes 
derivative with respect to time. 

We choose the scale factor such that aoHo = 2. If we 
are located at the origin, w — 0, then an object at redshift 
2 is at a comoving distance, 

dz' 



w{z) 



2 Jo E{z'y 



(38) 



and the comoving distance to the horizon (or the conformal 
time today) is 

dz' 



no = 7: 



2 70 



(39) 



4.2 Weak-Lensing Power Spectrum 

Given a power spectrum P{k, z) for the mass distribution, 
the power spectrum for the gravitational potential is 



P4,{k,z) = k~ 



'^{aoHofno{l + z)Yp(k,z) 



(40) 



In linear theory, the time evolution of the power spectrum 
is given by 

P{k,z) = P{k,z^ Q)[D{z)/D{z = 0)]^ 



where 



D{z) = 



5^0 E{z) 



l + z' 
[E{z')f 



dz' , 



(41) 



(42) 



is the linear-theory growth factor (see, e.g ., Peebles 19 93). 
The weak-lensing power spectrum is ( Kaiser 1992 ) 



P,,(^)=4 / J 



w{z) 



w(z) 



P^{n/w,z), (43) 



where g{z) is given in terms of the survey redshift distribu- 
tion, dN/dz, by 



w{z) 



dz 



I 1 dN 
N'dl' 



w{z) 



■w(z') 



Inserting Eqs. (|^-(^ into Eq. (^, we find 



where 



I4k)P{k,z = 0) dk, 



/vo 



In{k) = 144no[p(w = K/fc)]^[l + 2(w = /t/fc)]^ 

x[D{w = n/k)/Dof. 



(44) 



(45) 



(46) 



The mean-square variance in a cell with window function 
W{k) is 



where 



dkP{k,0)G{k), 



(47) 



G(fc) = 



fclJO 



^\w{n)fUk). 



(48) 



4.3 Model for the Spatial Density Power 
Spectrum 

For the power spectrum, we use 
2^2 



P{k) 



8 



-SH(k/2)"T^(kpMpc/hr), 



(49) 



where T{q) is the usual CDM transfer function, kp — kHo/2 
is the physical wave number, and F ~ f2o/i is given more 
accurately in terms of Q,oh and the baryon fraction fif, by 
Eqs. (D-28) and (E-12) in Hu fc Sugiyama (1 996). For the 
transfer function, we use ( Bardeen et al. 1986[ ), 

ln(l + 2.34g)/(2.345) 



T{q) 



[1 -f 3.89q + (16.1g)2 + (5.46(7)3 + (6.71g)4]i/4 



(50) 



If the po wer spectrum is norma lized to COBE, the am- 
plitude Sh is dBunn fc White 1997|) 



Snin, Qo) 



1.94X 10-"n-°-^«5-°-o^ 
X exp[a(n ~ 1) -I- b{n — 1)^]. 



(51) 



If primordial density perturbations are due to infiation, 
then there will also be a stochastic gravity-wave background 
which contributes to the COBE anisotropy with an ampli- 
tude dependent on the spectral index n. In this case, a — 1 
and b = 1.97. If we make no such assumption and suppose 
that the stochastic gravity-wave background is negligible, 
then a = —0.95 and b — —0.169. If we are uncertain of the 
gravity-wave contribution to COBE, then the COBE nor- 
malization above (with no gravity-wave background) will 
provide an upper limit to the true amplitude of the power 
spectrum. 

Alternatively, the power spectrum may be normalized 
at small distance scales through the cluster abundance which 
fixes as, the variance in the mass enclose d in spheres of 
radius 8 ft-i Mpc, to ag ~ (0.6±0.1)!r2o "'^^ ( |Vianafc Liddk 
199C). In terms of the power spectrum. 



^1=2^ / k^dkP{k) 



SjiikpR) 



(52) 



where R — 8h~^ Mpc, and ji{x) is a spherical Bessel func- 
tion. Since we are using ao 7^ 1, kp (rather than k) enters 
into the argument of the spherical Bessel function. 



4.4 Results for the FIRST Survey 

Fig. |l| shows the weight function (z) that enters into the 
calculation of the weak-lensing signal for the four FIRST 
redshift distributions considered by Cress & Kamionkowski 
(1997) and shown in Fig. 1 therein. Our best estimate for 
the FIRST redshift distribution is "DP7 (med-z)" derived 
from a radio-source luminosity function due to Dunlop and 
Peacock (1990), but we also include two other plausible es- 
timates from these authors, "DP7 (high-z)" and "DPI," 
as well as a redshift distribution derived from a luminos- 
ity function due to Condon (1984). For comparison, we 
also show the weight functions obtained from assuming all 
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0.012 



0.01 



0.008 



-0.006 ^ 



0.004 



0.002 




Condon 
.DP7 (med-7). 
.DP7 (high-7)^ 



redshift 

Figure 1. Redshift wcig ht functions g^{z) for the four FIRST 
redshift distributions considered in Cress & Kamionkowski 
(1997). Also shown are g'^(z) that would be obtained if all the 
sources were at redshift Zs = 1 or Zs = 2. 



1 — I I I I Mil 1 — I I I I Mil 1 — I II I MM 1 — I I I Mill 1 — I I I I II 




0.001 0.01 0.1 1 

(h-' Mpc-') 



100 



Figure 2. Three-dimensional COBE-normalized power spectra 
P{k) for the four models listed in Table 1. Models 1-4 are repre- 
sented by solid, dotted, dash, and dot-dash curves, respectively. 
Also shown (the light dashed curves) are the window functions 
G{k) needed for calculation of the mean-square ellipticity at 
0p = 1° and Bp = 6', and the window function (light dotted 
curve) needed for the calculation of erg. 



sources to be at a redshift of Zj, = 1 or Zs — 2. Below, we cal- 
culate the predicted signals with all four FIRST redshift dis- 
tributions to assess the uncertainty in the predictions from 
imprecise knowledge of the redshift distribution. 

Fig. ^ shows the spatial power spectra of the mass dis- 
tribution for the four models listed in Table |l|. The hght 
dashed curves are the (unnormalized) window functions 
G{k) needed to calculate the mean-square ellipticity for 
1° X 1° square pixels and 6' x 6' square pixels (obtained 
with our best estimate, DP7 (med-z;), for the FIRST red- 
shift distribution). The hght dotted curve is the window 
function needed to calculate as- These window functions 
illustrate that weak lensing probes power over a wide range 
of distance scales, and that the ellipticity in 1° x 1° pixels 
probes the power spectrum on larger scales than as- Here 
we used the DP7 (med-z) redshift distribution. The weak- 
lensing window functions would be shifted very slightly to 
larger scales if we had used the high-z DP7 redshift distri- 
bution and to slightly smaller scales if we had used the DPI 
or Condon redshift distributions. Figs. ^ and ^ show the 
(unsmoothed) weak-lensing power spectra and correlation 
functions for the four (COB& normalized) models listed in 
Table |l| (again, using the med-z DP7 redshift distribution). 

Table |^ lists the predicted rms gradient component of 
the ellipticity for several flat COBE-normalized and cluster- 
abundance-normalized cold-dark-matter models with and 
without a cosmological constant for 6' x 6' pixels, 30' x 30' 
pixels, and 1° x 1° pixels calculated with Eq. (^). We used 
the med-z DP7 redshift distribution for these calculations. 
For the window function, we use 



W{k) ^2Ji{x)/x with 



K9r 



(53) 



and 9p is the pixel size.|j In all cases, the rms ellipticity 
scales with the pixel size roug hly as {{jgf) oc with 
/3 = 0.3 — 0.4. We also list agOlo'^^ where erg is that obtained 
when the power spectrum is normalized to COBE. Note that 
if the power spectrum is normalized to asQ^'^"^ — 0.6, as in- 
dicated by the cluster abundance, then the rms ellipticities 
are 0.018 and 0.012 for 9p = 20' and 1°, independent of the 
model (except for the ilo = 1, h = 0.5 model which dif- 
fers negligibly). Weak lensing is due to perturbations in the 
gravitational potential (rather than the mass distribution), 
and the amplitude of gravitational-potential perturbations 
is fixed by the cluster abundance. This is why the weak- 
lensing power spectrum and correlations functions shown in 
Figs, g and ^ for the COBE-normalized Model 1 are so much 
higher than the others: this model predicts a value of crgf^o'^'' 
significantly larger than the others. 

We have checked that our calculations agree reasonably 
well with those of Jain & Seljak (1997) for the models and 
redshift distributions they consider. Their work illustrates 
that nonlinear effects (which we have not taken into account) 
are important only for Op < 10'. Corrections to the pre- 
dicted signal due to nonlinear evolution of the power spec- 
trum should increase the weak-lensing signal for 9p > 10', 
but only by a relatively small amount. 

Table □ lists the predictions for the COBE-normalized 
Qo = 1, h = 0.5, and n = 1 model for the four FIRST 



* Actually, this is the window function for circular pixels of the 
same area. The results should be similar if we use the window 
function for square pixels. 
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Table 1. Predicted RMS ellipticities, ((7g)^)^^^ of FIRST radio sources from weak lensing for Bp = 6', 20', and 1° for both COBR- 
normalizcd and cluster-abundance-normalized power spectra. We also list the values of eg for COB-E-normalized models. 



(100 ((7j)^)'^') 
COBE-normalized cluster normalized 



n 


h 


n 




0-8 




6' 


20' 


1° 


6' 


20' 


1° 


1 


0.5 


1 


0.0125 


1.21 


1.21 


5.2 


3.7 


2.2 


2.6 


1.8 


1.1 


1 


0.5 


0.8 


0.025 


0.71 


0.71 


2.8 


2.1 


1.4 


2.4 


1.8 


1.2 


0.4 


0.65 


1 


0.015 


1.07 


0.65 


2.6 


1.9 


1.3 


2.4 


1.8 


1.2 


1 


0.35 


1 


0.015 


0.74 


0.74 


2.9 


2.2 


1.5 


2.4 


1.8 


1.2 



Table 2. Predicted RMS ellipticities (in percent) of radio sources from weak lensing for the COBE-normalized Model 1 for 06' X 6' and 
1° X 1° pixels. 



Rcdshift Distribution 


100 ((7g^)2)'/' (ep = 6') 
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Figure 3. Weak-lensing power spectra Pgg(n) for the four 
(COBi5- normalized) models listed in Table 1. 



redshift distributions. If the high-z DP7 redshift distribution 
is adopted, rather than the med-z DP7 distribution, then 
the predicted weak-lensing signal is increased by about 6%. 
But if the true redshift distribution were more accurately 
represented by the DPI or Condon distribution, the signal 
would be smaller by about 20-25%. Note that the fractional 



10 ■='F 




9 (arcmin) 

Figure 4. Weak-lensing correlation functions C\(9) and C2{6) 
for the four (COBJ5-normalized) models listed in Table 1. 

uncertainty in the weak-lensing signal is smaller for larger 
pixel sizes than it is for smaller pixel sizes. 

4.5 Detectability of a Signal with FIRST 

Table ^ shows, for example, that for a COBE-normalized 
CDM power spectrum with flo — 1 and h = 0.5, the pre- 
dicted mean-square ellipticity in 1° square pixels is (0.022)^ 
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and (0.037)^ in 20' square pixels. The rms noise in 1° for the 
survey parameters used above {A = 10,000 deg^, n = 40 
deg"^ and e = 0.4) is (0.0075)^ and it is (0.015)^ for 20' 
square pixels, which gives signal-to-noise ratios of 9 and 6 
for 1° and 20' pixels, respectively. If the power spectrum is 
normalized to the cluster abundance, then the signal is just 
near the detection threshold. Since the signal increases with 
smaller pixel size as 6~^^ (with j3 = 0.2 — 0.3) and the noise 
increases as 9~^, the sensitivity decreases slightly with if a 
smaller smoothing scale is chosen. 

Ho wever, the sensitivity of the signal can be improved 



significantly if the full information encoded in the power 
spectrum is exploited with a maximum-likelihood analy- 
sis. For example, for the COBE-normalized f2o = 1 and 
h — 0.5 model, = 1.6 x 10"^" for the window func- 
tion corresponding to 20' pixels. From Eq. (|35|), the signal- 
to-noise with this maximum-likelihood technique would be 
A/<JA = 22, which is much larger than that obtained by 
just comparing the predicted and measured mean-square el- 
lipticity using either 20' or 1° pixels. Therefore, a proper 
maximum-likelihood analysis can improve the sensitivity by 
a factor of 2-3. Given that the signal for cluster-abundance- 
normalized power spectra is only on the verge of detectabil- 
ity when only the mean-square ellipticity is measured (i.e., 
the predicted mean-square ellipticity is only slightly larger 
than 3ct), we conclude that, with this more sophisticated 
maximum-likelihood analysis, a high-significance detection 
(> 6a) should be possible with these survey parameters. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have calculated the predicted ellipticity cor- 
relations of FIRST radio sources expected from weak gravi- 
tational lensing due to mass inhomogeneities along the line 
of sight for several plausible power spectra for the large-scale 
mass distribution in the Universe. We discussed the tensor 
Fourier analysis and statistical techniques needed to isolate 
the signal in the data. The shear field reconstructed from 
measured ellipticities can be decomposed into a "gradient" 
and "curl" component. Weak lensing predicts the presence 
of only a gradient component. Measurement of the curl com- 
ponent can be used to look for non-lensing artifacts in the 
data. 

We also estimated the amplitude of a signal which could 
be detectable with a survey as a function of the survey's 
source density, mean intrinsic source ellipticity, and area of 
the survey. We found that a detection of the signal from 
cluster-abundance-normalized power spectra could be ex- 
pected with good statistical significance (> 6cr) with sur- 
vey parameters which approximate those of FIRST. COB& 
normalized models produce an even larger signal. 

In addition to the statistical errors which we have taken 
into account, there will be systematic effects in the data 
which will mimic the effects of weak lensing. However, the 
most egreg ious of these effects can be modeled and cor - 
rected for ( [Refre gier & Brown 1997^; Refregier et al. 19971) 



result would place strict upper limits on the amplitude of 
mass (rather than luminous-matter) inhomogeneities in the 
Universe. Seljak (1997) has recently discussed application of 
more sophisticated statistical techniques developed primar- 
ily for clustering and the cosmic microwave background to 
weak lensing from large-scale structure. These will be needed 
for precise determination of the power spectrum for future 
weak-lensing surveys with better sensitivity. 

Although there are several other searches for weak- 
lensing correlations with optical surveys (e.g.. Mould et al. 



1994), as well as some recent claimed detections ( Villum- , 
sen 1995] ; Schneider et al. 1997), these optical surveys probe 



the ellipticity correlation function on much smaller angular 
scales than FIRST, which will probe the correlation func- 
tion on scales > 1°. Therefore, by combining the results of 
these surveys, the weak-lensing power spectrum can be re- 
constructed over a wide angular range. Although the signal 
may be more easily detected with optical surveys, these will 
probe scales where corrections due to nonlinear evolution of 
the power spectrum may be significant. On the other hand, 
FIRST will probe the power spectrum in a regime where 
nonlinear effects are small, so the comparison with theory 
will be less hampered by theoretical uncertainties from non- 
linear effects. 

The Sloan Digital Sky Survey (SDSS) will provide yet 
another data base with which to look for the effects of weak 



and it should be possible to approach the statistical limits 
discussed here. Even with a slight degradation of the signal- 
to-noise expected from systematic effects, the effects of weak 
lensing should be visible for cluster-abundance-normalized 
power spectra with a maximum-likelihood analysis. A null 



lensing on large angular scales (Stebbins et al. 1996). How- 
ever, SDSS sources will typically be at smaller redshifts. 
Therefore, by comparing results from the SDSS and FIRST, 
the redshift distribution of the weak-lensing distortions can 
be disentangled. Since FIRST and the SDSS will cover the 
same region of the sky, one can also cross-correlate the 
shear field indicated by FIRST with the foreground density 
field mapped by the SDSS. This will provide more stringent 
probes of the power spectrum and should also allow a direct 
measurement of the bias of SDSS sources. 
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